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In a perturbed Universe, comoving tracers on a two-dimensional surface of constant observed red- 
shift are at different proper time since the Big Bang. For tracers whose age is known independently, 
one can measure these perturbations of the proper time. Examples of such sources include cosmic 
events which only happen during a short period of cosmic history, as well as evolving standard 
candles and standard rulers. In this paper we derive a general gauge-invariant linear expression for 
this perturbation in terms of space-time perturbations. As an example, we show that the observed 
temperature perturbations of the cosmic microwave background (CMB) on large scales are exactly 
given by these proper time perturbations. Together with the six ruler perturbations derived in pQ, 
this completes the set of independent observables which can be measured with standard rulers and 
candles. 



I. INTRODUCTION 

Essentially all cosmological observations are based on 
detecting light emitted or absorbed from astronomical 
objects such as galaxies. From the direction that a pho- 
ton is observed and the shift in the frequency of the pho- 
ton, we infer the location of the emitter by assuming that 
the photon has traveled along a straight line; more pre- 
cisely, that it followed a geodesic of the homogeneous and 
isotropic background Universe. The actual path of the 
photon, however, is deflected from the straight line due 
to cosmic structures around the emitter as well as those 
along the line of sight from the emitter to observer. That 
is, the light path follows the geodesic in the perturbed 
universe, which is also perturbed from that in the back- 
ground universe. 

The deviation of photon paths from 'straight lines' 
leads to differences in the observed correlation functions 
of galaxies from the intrinsic ones. Recent studies [2r|5] 
have shown that the most dominant light-deflection effect 
comes from the scalar metric perturbations whose contri- 
bution to the galaxy clustering is negligibly small on the 
inter-galactic scales but induces a factor of few change 
on near-horizon scales. In particular, this effect shows 
the same scaling as the scale-dependent bias signature 
due to primordial non-Gaussianity and correspondig to 
A/nl < 1; thus it has to be correctly modeled for future 
galaxy surveys which is pursuing the non-Gaussianity pa- 
rameters with similar accuracy |S] . 

On the other hand, the same light-deflection can also 
be used for studying the clustering and growth of cos- 
mic structures. The most popular method along this line 
for studying large scales structure is weak gravitational 
lensing. Here, the primary observable is the coherent 
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structure, or clustering, in the cllipticity of galaxies on 
large scales. Because intrinsic correlation of ellipticities 
of galaxies on large scales is expected to be very tiny, we 
can attribute the measured correlation to the correlation 
in the light-deflection due to cosmic structure. 

Beyond the conventional weak lensing method, p] have 
shown that observables such as the length of standard 
rulers and the luminosity of standard candles are system- 
atically distorted by the light-deflection and can thus be 
used as proxies for the large-scale cosmic structure. The 
distortion of an (intrinsically) spherical object has six in- 
dependent components which are scalar (2), vector (2), 
and tensor (2) under the rotation on the celestrial sphere. 
Note that one of the scalar modes and two tensor modes 
are the standard weak lensing observables: magnification 
and shear, respectively. In pQ, we have presented a co- 
variant formalism for these six components in terms of 
the metric perturbations and peculiar velocities. 

In this paper, we shall study yet another observable 
that is distorted by the light-deflection: cosmic clocks. A 
cosmic clock refers to a spacetime event with observable 
proper time since the Big Bang, as measured by a comov- 
ing observer. That is, any global event with which we can 
synchronize a space-like hypersuface in terms of proper 
time is a candidate cosmic clock. The examples of the 
cosmic clock includes BBN (Big-Bang Nucleosynthesis), 
last scattering of cosmic microwave background (CMB), 
thermal decoupling of CMB photon and baryons, begin- 
ning and end of the reionization, etc. Another class of 
cosmic clocks can be set by using time evolution of ob- 
servables such as the mean number density of a certain 
type of galaxies, or the length of a time-evolving cosmic 
ruler such as the physical size of galaxies. In this case, 
having an observable proxy for the proper time, we can 
reconstruct the hypersufaces of constant proper time. 

Each cosmic clock event may be identified through 
various observational signatures, but in general this will 
involve detecting light from some source. Most impor- 
tantly, then, we can measure the redshift of the photon 
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from each cosmic clock event. Because of the pertur- 
bation to the photon geodesies, however, the measured 
rcdshift of photons emitted cosmic clock events at a fixed 
proper time varies over the celestial sphere, and it is this 
variation that we shall study in this paper. We can of 
course cquivalently phrase the variation is redshift per- 
turbation from a constant-proper-time slice, or proper 
time perturbation on a constant-observed-redshift slice. 

This paper is organized as follows. After deriving the 
gauge-invariant formalism for the proper-time perturba- 
tion in § [nj we present two examples where the proper- 
time perturbations become important: the scalar- type 

and super- 
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distortions for evolving standard rulers (§ [III ) 
horizon temperature anisotropies of the CMB (§ IV). Wc 
conclude in §[V]with discussion. App.[X]proves the gauge 
invariance of the proper-time perturbation, and App. [C] 
contains explicit expression of the proper-time pertur- 
bation in terms of the density contrast in synchronous- 
comoving gauge, useful for performing quantitative cal- 
culations. 



II. FORMALISM 
A. Notation 

We write down the most general form of perturbed 
FRW (Friedmann-Robertson- Walker) metric as 



A|| = n t X\ 
E\\ = fiifijE 13 , 

x]_ = r ij Xj 



(5) 



Correspondingly, we define projected derivative opera- 
tors, 



<9|| = n l di, and 
<r Vi>,. 



(6) 



Note that d]_, dn and d±, do not commute. Further, 
we have 



d j n i =d ±j n i = -V, 1 , 



(7) 



where x is the norm of the position vector so that n l — 
x l /x- Note that n l and dn commute. More expressions 
can be found in § II of [5J. 

For all numerical results, we shall assume a flat ACDM 
cosmology with h — 0.72, O m = 0.28, a scalar spectral 
index n s = 0.958 and power spectrum normalization at 
z = of erg = 0.8, which is consistent with cosmological 
parameters estimated from WMAP9 [S], and reasonably 
close to the results from Planck [TO] . 



ds 2 = g ^ u dx^ dx u 

= a 2 (r/) [ - (1 + 2A)drf - 2B i di 1 dx i 

+ (Sij + hij) dx l dx 3 



(1) 



where we have assumed that the background Universe is 
spatially flat. Here, r\ denotes conformal time and a{rf) is 
the scale factor. Following usual convention, the spatial 
part is further expanded as 



2D8 ij +2E ij , 



(2) 



where Eij is a traceless 3x3 tensor. We shall also present 
the end results in two popular gauges: the synchronous- 
comoving (sc) gauge, where A = = Bi, so that 



ds 2 = a 2 (rj) [-drj 2 + (5y + hy) dx l dx j ] 



(3) 



and the conformal-Newtonian (cN) gauge, where Bi = 
= E^. In the latter case, we denote A = "J, D = $, 
conforming with standard notation, so that 

ds 2 = a 2 {ri) [-(1 + 2^)drf + (1 + 2<&)6 ij dx i da?] . (4) 

It is useful to define projection operators parallel and 
perpendicular to the observed line-of-sight direction n l , 



B. Proper time perturbation 

Consider the redshift perturbation of the set of cosmic 
clock events defined by a constant proper time tp. We de- 
fine the "cosmic clock" observable T(n) as the difference 
in In a between a constant-proper-time surface tp = const 
and a constant-observed-redshift surface z = const. Al- 
though phrased as a perturbation in In a, we will fre- 
quently refer to T loosely as the proper time perturba- 
tion. This is because at leading order, the perturbation 
to the proper time Atp(n) at observed redshift z is sim- 
ply related to T through 



At F (h) = H-\z)T{h) 



(8) 



Note that, since it is defined by two observationally 
well-defined quantities (proper time and observed red- 
shift), the perturbation T is clearly an observable; thus, 
whatever expression is obtained for T has to be gauge- 
invariant. 

The proper time interval dtp is defined through 



dtp = y / — g^dx^dx" . 



(9) 



A comoving source with velocity v % obeys (in comoving 
coordinates, Eq. and at linear order in v) 



dx a = d-q; dx l = -v l dt = v l dr\ . (10) 
a 
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We then have, to linear order in perturbation, 



expression for the perturbation T: 



dtp = (1 + A)adrj. 



(11) 



Integrating Eq. (11), we obtain an expression for iF|?j,x> 
the proper time of a comoving source passing through x 
at coordinate time rj, at linear order 



£f U,x — 



[1 + A(x,rf)]atf)drf. 



(12) 



In the case at hand, rj is the coordinate time of emission 
of the observed photon. The ratio of scale factors at co- 
ordinate time r\ and at the proper time tp of an observer 
passing through (77, x) is then given by 

= 1 + H(rj) I A(x, rj')a(r)')dr)' . (13) 



Here, a[fj(tp\ 7hx )] denotes the scale factor in an unper- 
turbed Universe at the proper time tp that a comov- 
ing source has when passing through the spacetime point 
fox). 

Let us consider a standard ruler whose proper length 
evolves in time. Then, by using Eq. (13 1, we can 



parametrize a time evolution of the proper size of the 
standard ruler ro(a) through its value in an unperturbed 
Universe as function of the scale factor a. We can then 
write the ratio between the actual proper size of the ruler 
ro(a(tp\ x o x )) and the size of the ruler if the proper time 
of emission coincided with the age of the background Uni- 
verse corresponding to the observed redshift z, 



f"o(a(*FU°, x )) _ r (a (t F \ x a^)) r {a{x )) 



ro(a) 



r (a(x )) r (5) 



= 1 



din ro(d) 
din a 

din ro(d) 
din a 



In 



(*fU°,x) 

a(x°) 



In 



a(x°) 



T. 



(14) 



Here, x° is the coordinate time at which th e ph oton was 
emitted, and a = (1 + z)^ 1 . Note that Eq. (14) assumes 
that a a = 1 at observation (z — =>■ d = 1), i.e. r (l) 
corresponds to the ruler scale today as calibrated by the 
observer. This clearly implies that T — for a locally 
measured ruler (any non-zero value would be merely a 
constant offset and could be absorbed into r ). Note 
that the epoch of observation to is really fixed in terms 
of proper time, rather than coordinate time. 

We now use the fact that \n[a(x°)/a] is precisely the 
perturbation A In a derived in [T] . Note that A In a is not 
gauge- invariant itself. With this, we arrive at the explicit 
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T = H A[x,r]']a(r)')dr) ' + Alna 
Jo 

rn rVo 

= H A[x,rf]a(ri f )drf - H A[0,t]']a(r]')drf 
Jo Jo 

+ A - A + v\\ -v\\ + [ dx 
Jo 



(15) 



A + l -h\ +B' n 



where 77 is defined through a(fj) = a, and x = Vo ~ V 4 l n 
the second line, quantities without a subscript are eval- 
uated at the source, while quantities with a subscript o 
are evaluated at the observer. The terms under the x 
integral are to be evaluated along the photon geodesic. 
Further, H — H(d), and we have used the expression for 
A ln a derived in £Q : 



Alna = A n - A + vu - vu 



x 



dx 







- A ' + l h \\ +B \\ 



H / A(O t r/)atf)dri 
Jo 



(16) 



The last term ensures that the observer resides at a fixed 
proper time. This term is the only non-vanishing contri- 
bution to A ln a in the limit z — > 0, so that Eq. ( 15 1 yields 
T —> in this limit as desired. T has two sources: the 
perturbation of the apparent coordinate time of emission 
due to Dopplcr shift, gravitational redshift and ISW (In- 
tegrated Sachs- Wolfe) effect, and the perturbation to the 
coordinate time at fixed proper time of the source. 

We again emphasize the difference between T in 
Eq. (15) and Alna in Eq. (16): while Alna gives the 
perturbation of a z — const surface from a constant- 
coordinate time surface, the perturbation T is the pertur- 
bation of the z = const surface from a constant-proper- 
time surface. Since both z and the proper time of a source 
since the beginning of the Universe are observable, T is 
observable, while Alna is not. The gauge-invariance of 
T is shown explicitly in App. \K\ for scalar perturbations. 

We can construct an explicit procedure for observing 
T(n) as follows. Consider two rulers a, b which are scaled 
so that r 0a — r fc at some proper time tp but which evolve 
differently (dlnro a /dhia ^ dlnrob/dlna). Assume fur- 
ther that when averaged over the sky (or the survey area) , 
this proper time corresponds on average to a redshift z. 
In that case, for sources at redshift z, all projection ef- 
fects drop out in the local difference r$ a — measured in 
a given direction n. The only remaining contribution to 
the apparent difference between the rulers a and b then 
is 



rpg - r ab 
roa 



d In r 0a d In r b 



din a 



dim 



T(n). (17) 



In general, T is one out of several contributions to the 
observed perturbations of standard rulers, as discussed 
in § |III| Another example is a set of sources which emit 
at a fixed proper time tp. Then, —T corresponds to the 
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perturbation in observed redshift of these sources. We 
will consider this case in § |IV| Note that, as for the other 
ruler perturbations (see App. C in pQ), the lowest order 
contribution of a Fourier mode with wavenumber k is of 
order (k/ Hq) 2 . This is shown in App. pi 



For convenience, we give the expressions of Eq. ( 15 ) 
specialized to synchronous-comoving gauge, 



1 

2 Jo 



h '\\(nx,Vo-x) dr]' , 



(18) 



and conformal-Newtonian gauge, 

(T)cN =H *[x, v'Hv'W - H / * [0, rf]a{rf)drf 
Jo Jo 

rx 

+ * - \P + V\\ - V\\ 



[ d X [&-V] A> 

Jo 



(19) 



Note that (T) sc = (Alna) sc = 5z, with Sz as defined 
in This is because in synchronous-comoving gauge 
the coordinate time coincides with the proper time of 
comoving observers. 
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III. EVOLVING STANDARD RULER 



Using the results of the previous section, it is now 
straightforward to generalize the case of a fixed standard 
ruler considered in [T] to an evolving ruler. Specifically, 
Eq. (30) in that paper becomes 



rl{a) 



r 2 



2T 



d In ro (a) _ 
din a 



2 + 2Alnar 2 



+ 2a 2 ( v\\ Sx 2 + vj_ iSx l j_5x\\ J 

+ 2& 2 5 i:j Sx t (5x^8% + Sx^d^k) Ax j . (20) 

Here, f is the apparent size of the ruler, while ro(a) 
is the true size of the ruler in an unperturbed Uni- 
verse evaluated at the apparent scale factor at emission 
a = (1 + z)^ 1 . Thus, given knowledge of the ruler as 
function of time ro(a), we can measure the individual 
contributions to Eq. pO] ). Moreover, the evolving ruler 
case is probably more common than a fixed ruler, when 
applied to sizes of galaxies, correlation lengths of a tracer, 
or the BAO feature (which is fixed in comoving coordi- 
nates, ro(a) cx a). 

The additional term cx T does not spoil the decompo- 
sition of [I] into parallel and perpendicular components 
relative to the line of sight, 



ro 



-C 



(Sx») 2 



■Bi 



5x\\Sx\ 



, Sx\ 8x\ , . 
Aij^^, (21) 



where f c = r/a is the apparent comoving size of the ruler. 
Rather, using Eq. (20) and Eq. (21) we can immediately 



FIG. 1: Angular power spectrum of T (blue solid) for scalar 
perturbations in the standard ACDM cosmology (§ II A I. For 
comparison, we also show the power spectra for the magni- 
fication M (magenta dash-dotted) and longitudinal scalar C 
(black dashed), calculated for a non-evolving ruler. All quan- 
tities are evaluated for a fixed source redshift of z = 2. 



read off the contribution to the longitudinal component 
C, 



<iln?'n(a)^ 1, „ . , 

C = ^^T- Alna h n -v\\ - d-sAxu (22) 

din a 2 11 11 x 11 



and to the magnification M. , defined as the trace of A4 



M = V ij Ai 



:2 dlnro(a) 1 ( , 

dlnd 2 V 



2k 



-Aim 



X 



(23) 



As a scalar on the celestial sphere, T does not contribute 
to the vector Bi and the transverse components of .Ay 
(shear) . 

As derived in App. [C] Cj-(£) is given in terms of the 
matter power spectrum today P m {k) (in synchronous- 
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comoving gauge) by 



C T {1) = - / fc 2 dfcP m (fc)|F/(fc)| 



(24) 



F< (k) = H 



([g-l]D„_a) kr dr ]]l {x) + F^ a {k) 



F, Alna (£;) 



aHJD 1. . ... 
— ^— d £ --(g-l)D*_ ) ji(x) 



dxDiswji(x) 



Here D is the matter growth factor and £>$_ oc fc~ 2 
is the relation between matter and potential perturba- 
tions in cN gauge (see App.[C]for details). The quantita- 
tive importance of T is illustrated in Fig. [T] which shows 
the angular power spectrum Cf(i) of T at fixed redshift 
5 = 2. We see that T is significantly smaller than A4 
and C except on the very largest scales. The reason is 
that the contributions to T are suppressed with respect 
to the leading contributions to Ai and C by a factor of 
aH/k, where k is the typical wavenumber contributing 
to a given angular scale. In fact, for this (and smaller) 
source redshift, Cf(l) is completely dominated by the 
peculiar velocity contribution. The typical wavenumbers 
contributing to A4 at I ~ 500 are at the turnover scale 
of the matter power spectrum, k ~ 0.01 h/Mpc. We thus 
expect that Cj-(£) is suppressed with respect to Cm{^) 
by (aH/k) 2 ~ 10 3 at those scales, which roughly matches 
the numerical result. 

While the exact relative contribution depends some- 
what on the source redshift, the order-of-magnitude sup- 
pression remains the same. This is illustrated in Fig. [2] 
which shows Cj-(Z) for a wide range of source redshifts 
up to z — 100 (at / < 100, the gree n dashed line in fact 
shows Cf(l) for z ~ 1100, see § IV). The dotted lines in- 
dicate the angular power spectrum of the matter density 
contrast 5^ on constant-proper-time slices. This serves 
as an illustration of the typical amplitude of intrinsic 
perturbations to tracers. Clearly, T is subdominant ev- 
erywhere apart from the largest scales at high redshifts. 
Note however that the power spectrum from 5!£ can be 
suppressed if projected over a broad redshift range, as 
illustrated by the magnification in Fig. [I] 

Thus, unless one is dealing with a very rapidly evolving 
ruler (<91nro/<91na 3> 1), the contribution to the magni- 
fication and longitudinal ruler perturbation provided by 
T will be very small numerically. However, we next con- 
sider an example where the perturbation T turns out to 
be the dominant observed contribution on large scales. 



IV. CMB IN THE SUPERHORIZON LIMIT 

The observed CMB photons originate from the last 
scattering surface, which occurred at a fixed physical 
age t* of the Universe, that is, at constant proper time 
tp = t* for comoving observers. The value of is 
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FIG. 2: Angular power spectrum of T(n) (blue solid) and 
the matter density perturbation <5„ in synchronous-comoving 
gauge for different sharp source redshifts. From top to bot- 
tom, the curves show z = 1, 2, 5, 10, and 100, respectively. 
5^ serves as a rough order-of-magnitude estimate of intrin- 
sic tracer density perturbations on a constant-proper-time 
slice. The dashed green line near the bottom shows the power 
spectrum of the fractional CMB temperature perturbation & 
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obtained by combining atomic physics with the mean 
observed temperature of the CMB today (along with 
some assumptions about the stress-energy budget of the 
early Universe). The CMB temperature perturbations on 
scales that were superhorizon at recombination (i < 100) 
originate entirely from projection effects; in other words, 
the large-scale CMB temperature perturbations can be 
seen as a special case of the evolving ruler described 
above. Essentially, the standard ruler is in this case given 
by the photon occupation number I^/is 3 . 

Since we are dealing with a projected quantity, only 
the transverse perturbations Aij are relevant. Moreover, 
the CMB temperature is a scalar variable, so that the 
trace-free component of Aij does not contribute. We are 
left with the magnification M.. However, all effects of the 
propagation of light leave the photon phasespace density 
Ivjv 3 invariant (surface brightness is conserved), so that 
the non-evolving-ruler magnification of pQ does not con- 
tribute. Thus, the only contribution to the fractional 
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temperature perturbation O(n) is given by the proper 
time perturbation: 



T 



din a 



(25) 



where have used that T oc a 1 for a free-streaming black- 
body HQ. 

We now show this more explicitly. By definition of 
A In a, the relation between scale factor at emission and 
observed redshift z is given by 



l (0=(l + 5)- 1 (l + Alna). 



(26) 



The CMB temperature T(n) and observed redshift z are 
related by 

^ = [l + e(n)]^ = (l + 2)- 1 , (27) 
Jo Jo 

where T is the mean observed CMB temperature, To is 
the temperature at emission (essentially set by atomic 
physics), and O(n) is the temperature perturbation 
which wc intend to derive. The coordinate time at emis- 
sion x° is set by the requirement that it correspond to a 



fixed proper time tp = t*. Applying Eq. ( 13 1 to this case 
yields 



ft* 



1 + 



A[xL,rj\a(r))dr), 



(28) 



where 77* is the conformal time corresponding to t* in the 
background, a* = 0(77*), and £/* = H(r)*). Inserting this 
into Eq. ( 26 1 and using Eq. ( 27 ) , we obtain 



A(x, 77) adrj 



[l + e(n)]^(l + Aln a ), 
(29) 



which can be solved for to yield 
6(n) = - H* [ A(x,ri)a(r))dr)-Alna=-T(ti) 







(30) 



where we have used that T = a* To by definition of 
a* fa 1/1089. This expression is identical to — T [Eq. ( 15 ) 
for a = ft*l, thus validating our considerations lead- 



ing to Eq. (25). Up to a sign, T is the general, lin- 



ear, gauge-invariant expression for the CMB tempera- 
ture perturbation in the superhorizon limit, i.e. with- 
out acoustic contributions. In particular, adopting the 
conformal-Newtonian gauge [Eq. Q], and using the su- 
perhorizon limit where $(x) ps const, Eq. (30) reduces to 
[see Eq. (ph] 



0(n) 



1 



d x W - *'] 



H / *(0,77 )ft(ry)d77 
Jo 



(31) 



Here we have also used that recombination happened 
long after matter-radiation equality so that H„t„ =2/3. 
This is the well-known expression for the large-scale 
CMB temperature perturbation in conformal-Newtonian 
gauge. The terms on the second line only contribute 
to the monopole and dipole of 0(n). Thus, the large- 
scale CMB temperature perturbations are nothing else 
than minus the proper-time perturbations on a constant- 
redshift surface. Of course, at second order lensing deflec- 
tions do modify the statistics of the CMB temperature, 
an effect which can again be addressed in this formal- 
ism by using the intrinsic (Fermi frame) CMB correlation 
function as ruler |12) . 

Note that all these projection effects are independent 
of the photon polarization. They thus do not affect or 
induce polarization in the CMB. Instead, the polariza- 
tion is imprinted by the physical effects of the long- 
wavelength perturbation, and is correspondingly sup- 
pressed by (fc/ft*7J*) 2 in the low- A; limit, while the pro- 
jection effects T which determine the CMB temperature 
scale as (k/Ho) 2 in that limit (App.pl). 



V. DISCUSSION 

In this paper, we present the gauge-invariant ex- 
pression for the proper time perturbation on a two- 
dimensional surfaces of constant observed redshift. The 
proper-time perturbation can be measured from observ- 
ables which define the constant proper time hypersurface. 
We present two observables which allow for a measure- 
ment of this perturbation: the large-angle temperature 
perturbations of the cosmic microwave background and 
standard rulers with evolving proper length. More gener- 
ally, these two examples represent two classes of proper- 
time observables. 

One class of observables consists of cosmic events de- 
fined by a unique time and sufficiently short duration 
so that the proper time is well-defined. This class of 
events includes for example the epoch of neutrino decou- 
pling, Big-Bang Nucleosynthesis, CMB last scattering, 
thermal decoupling of baryons from the CMB, beginning 
and end of reionization, etc. As shown for the case of the 
last scattering surface of CMB, however, the proper-time 
perturbation we have calculated here dominates only on 
superhorizon scales at the epoch of emission. On smaller 
scales, other physical effects generate perturbations in 
associated observables. Because these small-scale per- 
turbations are typically of order of the density contrast, 
while the contribution from the proper-time perturbation 
involve velocities and potentials, the latter are relatively 
suppressed compared to the former on subhorizon scales 
as shown in Fig. [2j The angular anisotropics in the cos- 
mic neutrino background, which was emitted at the neu- 
trino decoupling epoch, are also given by the proper-time 
perturbation on scales that were superhorizon at that 
time. 

Another class consists of observables with known time 
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evolution which allow us (in principle) to define a con- 
stant proper-time hypersuface. This class of events in- 
cludes time evolution of standard rulers, time evolution 
of mean number density of specific population of galax- 
ies, etc. We have shown that the time evolution of the 
standard ruler alters the expression for radial distortion 
C and magnification M.. The amplitude of proper-time 
perturbation, however, is about three and one orders of 
magnitude smaller compared to C and Ai, respectively. 
Therefore, the proper-time perturbation is probably only 
important for rapidly evolving standard rulers which can 
make up for the difference in amplitude. 

The time evolution of the galaxy number density af- 
fects the observed clustering of galaxies [5J. The proper- 
time perturbation T here in synchronous-comoving gauge 
is equal to Sz in [BJ , and the time evolution of the galaxy 
number density yields a contribution to the observed 
galaxy overdensity of 



Appendix A: Gauge-invariance of T 



S„ D b e T ■ 



(1 + 5) 



<iln(a 3 n g ) 



dz 



T. 



(32) 



Again, this effect is suppressed compared to contributions 
such as the intrinsic galaxy density contrast, redshift- 
space distortion and magnification bias by k/aH, because 
T is dominated by the line-of-sight peculiar velocity. On 
the other hand, using multiple populations of galaxies 
[13) may help measure the proper-time perturbation by 
suitable optimal weights assigned to tracer densities in 
a given volume. Note that diffuse backgrounds of any 
wavelength can also be useful for this purpose since they 
are not affected by lensing bias. Multiple "tracers" can 
be implemented for example by thresholding. 

Finally, we point out that the proper time perturba- 
tions derived here can be seen as a test of homogeneity of 
the Universe [THUS] ■ That is, a measurement of, or upper 
limit on, the magnitude of 7~(n) consistent with the nu- 
merical results for ACDM presented here would provide 
direct evidence for the assumption that the unperturbed 
background FRW metric provides a good description of 
the observed Universe, i.e. that the Copernican principle 
holds. 



We now confirm the gauge-invariance of Eq. (151 under 



a general scalar gauge transformation. Writing 



d l L{x^) 



(Al) 



we use the transformations of the metric perturbations 
given in App. Al of [6j. Relevant transformations for the 
case at hand are 

A -)• A - aHT - T 
v ->• v + L' 
B^B + L' -T 
<p —tip — aHT 

7 ->7 — L , (A2) 

where we consider only scalar modes of the metric per- 
turbations 



Bt = B, 



2ip6ij + 27,., 



(A3) 



In terms of the scalar perturbations, the line-of-sight pro- 
jection of the metric perturbations are 



B\\ = n l B^ 



(A4) 



Under gauge transformation, 



1 



- A + -h\\ + B\\ = -A + tp + n l n J 7 j y 



n l B 



A + aHT + T' + ip— aHT + n i n j % ij - u'n'l.,, 
h i B i + fi b L', 



A + ip + h l n 3 ^ A j + h % B A 
T — n l n J ' L t . 
1 



n'i „■ — n l T 



= - A + -h\\ +B\\ -d x (T + d\\L) 
Here, we have used that d x — d\\ — d v , and 



(A5) 
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= 9||(L' - d {l L) = -d x d l{ L. 
Using above, we find that 



(A6) 
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dx 



-A' + -h\ + B[ 



d X 



-A 1 



1 



h'n + B'u 



(A7) 



-T'{x)-d\\L'{x) + T' + d\\L' . 

Combining this with 

A Q - A + U|| - U|| -> A Q - A + U|| - 
- a H T - T + aHT + T + a, L' - a, L' a (A8) 



yields 



A + vn 



->A„ — A + vu 



x 

"X 

'•ii.-. / dx 





-A' + -h\ +B[ 



-A' + -h\ +B[ 



(A9) 



Using Eq. ( 15 ), we obtain the gauge transformation of T 
as follows: 

T(n) -> T(n) + &HT[xh; fj] - a o H o T[0; Vo ] 

-H {aHT[xm;fj(t)]+T'[xm;v(t)]}dt 

+ H o J " [aHT[0]ij(t)]+T , [0;fj(t)]}dt 

= T(n) + aHT cm - a H T - H (aT em - [aT] a=Q ) 

+ H a (a T - [aT] o=0 ) 
= T(A). (A10) 

Here, we have abbreviated T em = T[xn; 77], T G = T[0; ?/ ], 
and have further assumed that T remains finite as a — > 
(or at least diverges less rapidly than a -1 ), so that we 
can neglect (aT) evaluated at a = 0. Thus, T defined 
through Eq. (15) is gauge- invariant as required for an 



actual observable. 



Appendix B: Pure-gradient metric perturbation 

It is instructive to show that a constant+pure gradi- 
ent metric perturbation does not contribute to T. In 
this case, the Sachs- Wolfe term cancels exactly with the 
Doppler redshift, and the lowest order contribution by a 
single Fourier mode of wavenumber k is proportional to 
(k/H ) 2 , as pointed out by [TCI IT?]. 

For simplicity, we specialize to an Einstein-de Sitter 
(EdS) Universe where distance and growth calculations 
are particularly simple. In EdS, the linear growth factor 
becomes D(a) = a, and we have $ = — V?' = and 



Ht = 2/3. In this case, Eq. (19) becomes 



(T)cN = S \ [* - *o] + *o - * + «|| - «|| 

= i[*(0)-*(xn)]+«||-«|, o . (Bl) 

We now consider a constant+pure gradient potential per- 
turbation, 



9(x,rj) = ¥ [1 + k-x] 



(B2) 



where \&o and k are constants. As before, the observer 
is assumed to be at x = and to be comoving. We then 
obtain 



1/2 T 2 h T 

v = a 1 ' 2 — V\\ n = L * - 



(B3) 



The EdS background yields 

1,1 da 2 



X : 



g a 2 H(a) H 



(B4) 



The monopole Q(k ) contribution to T is clearly vanish- 
ing from Eq. (Bl). The dipole component oc k\\ is given 
by 



T' 



^ fc HX*o + ^(l-a 1/2 )*o = 0, (B5) 



as desired. Note that the observer terms in Eq. ( 19 ) 



are crucial for obtaining this result. The lowest order 
contributions then appear when expanding Eq. (B2) to 
quadratic order in k, and scales as (k/H ) 2 - 



Appendix C: T in terms of synchronous-comoving 
matter density perturbation 

In this section we derive T in terms of the familiar mat- 
ter density contrast 5^ in synchronous-comoving gauge. 
We make use of the relations in App. F of [1 . In partic- 
ular, we write the potential ^ in Fourier space as 

*(k )77 ) = ^[ 5 -lp $ _^(k, ?7 o). (CI) 

In a ACDM cosmology (or more generally for a smooth 
dark energy component), we have 



-£>*_ (k, 77) = 3n m ^^D(a(r))) 



a?H 2 
Hi 



:3a„ -^a- 1 (7 ? ) J D(a(7 7 )) 



5(M)=0, 



(C2) 



where D(a) is the matter growth factor normalized to 
unity at a = 1. Here, a subscript denotes that the 
quantity is defined at the present epoch rj = 770, while 
a tilde denotes quantities evaluated at the inferred scale 
factor at emission a = (l + z) -1 . We will denote the 
power spectrum of at z = as P m (k), and define 
x = kx-i x = kx- Using that 

T = H [%[x,r)']a(ri')dri' + (Alna) cN , (C3) 
Jo 

the contribution of a single Fourier mode with wavevector 
k to T(n) is given by 



T(k,n) =6 



1 



([g-l}D^_a) kv dr 1 e^S^(k,7 1 o) 



(Alna) cN (k, n) . 



(C4) 



where fi = k • n and (Alna) C N(k, n) was derived in 
App. Fl of [JJ. The angular power spectrum of T(n) 
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(at fixed observed redshift z) is then given by 



We now use the equation for the growth factor D(rf), 



Cr{l) = -J k 2 dkP m (k)\Ff(k)\ 2 
Ff(k) = h£ \ {[g - 1]D»_ a) ^ dr, ji{x) 



F l Alna {k) 



(C5) 



aHfD a 1 , 



(C6) 



which yields (aD 1 )' = m m0 H' Q 2 D/2. Thus, 



(Cll) 



where 



dz Di SW (k,z)ji(x) 







(C7) 



Asw(M) = ^^.(MO))- (C8) 

Here we have neglected pure monopole and dipole con- 
tributions which are straightforward to include but not 
relevant observationally for all measurements considered 
in this paper. 



We now restrict Eqs. ([C6|)-(|C7| to ACDM. Using 



Eq. (C2|, we have 



' d X 3H 3 a 3 n m (a)D(a) [f(a) - l]ji(k X ) 

(C9) 



where f(a) = d In D(a)/d In a. The first term in Eq. (C6) 
on the other hand becomes 



6 [\{\9-i\ D *-*) k « 



H 



s _ n d,T]ji(x) 

n m0 ^D 3l (z). (Cfo) 



H l V ('I) ^H 2 Q D = -H f{ajy)' = -H[aD')l 



-(aB)< 



dD 



e?lna 



-(dH) 2 fCa)D(fj). 



(C12) 



Thus, the time integral term acquires the surprisingly 
simple form 



and we obtain 



k 2 



-fDji{x), 



(C13) 



pfik) AC = DM 



(C14) 



+ / rf X 3F 3 a 3 fi m (a)D(a)[/(a)-l]jKfc X ). 
Jo 
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